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CHAPTER I 


INTRODUCTION 


It has long been recognized by meteorologists that there 
is an interdependence between the pressure distribution and 
wind flow on a level surface. In the steady state the pres- 
sure, Coriolis and centripetal accelerations all sum vector- 
ially to zero. The steady state is a condition which may be 
frequently approximated in nature, but it is probably seldom 
realized. Thus, there is a process of continual adjustment 
going on in the atmosphere. The wind accelerates in response 
to changes of pressure gradient, and pressure gradient changes 


in response to non-gradient winds. 


Rossby [ 8] and Cahn [3] studied the mutual adjustment of 
pressure and velocity distributions in certain simple current 
systems. These currents were narrow but infinite in length, 


having a current profile similar to that which has now become 


associated with the jet stream. 


Unbalanced forces appear to occur in the atmosphere in 
another important fashion, which has been emphasized recently 
by O'Connor [6 | because of its importance in prognostic anal- 


ysis. In this situation strong winds are found either upstream 





or downstream of light winds. This condition is observed 

daily between the isotach maxima and minima occurring along 

the jet stream on the upper-level facsimile charts prepared by 
the WBAN Analysis Center. Under these circumstances it is easy 
to see that even though all forces may be balanced at a given 
instant, inertia may carry air into a region with speeds which 
are not compatible with the pressure distribution existing 
there. What occurs then is not well understood, but there 

must be some mutual adjustment between the wind velocity and 
the pressure distribution. 

Brunt i) discusses an investigation of this problem by 
Brunt and Douglas. This treatment of the problem makes a nun- 
ber of simplifications and does not show the effect of the 
curvature of the air motion. Fetterssen [7] discusses this 
problem te a slightly different manner, but makes the same 
approximations as Brunt and Douglas. The method of Petterssen 
is followed and extended in this study, but the equations of 
motion will be expressed in terms of the topography of an iso- 
baric surface instead of the pressure distribution ina level 
surface. 

All equations and the discussion which follow refer spec-— 


ifically to the Northern Hemisphere. 





CHAPTER II 
DERIVATION OF THE ISALLOHYPTIC GRADIENT EQUATION 


The equations of frictionless horizontal flow on an iso- 


baric surface referred to a natural coordinate system are 
: £ 4 SE 


~ 
and ay = -a<F (2.1) 


These may be written in vector form as follows: 


ey) Sey 
g KX, z =(f{tkv) VE In (2.2) 


Differentiating partially with respect to time, we obtain 
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Taking a cross product of both sides with H gives 
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This is a form of the Brunt Douglas equation. The left 
side, excevt for the factor pet is the negative of the isallo- 
hyotic gradient. No terms have been ignored, and the only 
assumption thus far is that the motion is horizontal and fric- 
tionless. 

In this form (2.3) is difficult to evaluate, since every 
term on the right-hand side involves a time derivative. Other 
expressions may be obtained which may be substituted for some 
of these. One of these is obvious from the tangential component 
Prot 2.1). Thus -0-2 may be substituted for wy in the first 
term on the right of (2.3). 

This type of direct substitution is not valid in the term 


3 


/ 
ae however, since the coordinate system used is not fixed 


(see Figure 1). In this case 
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where /’ is the angle between the / axis of a fixed // coordinate 
system and the unit vector 7 which is directed along the 
positive © axis of the natural coordinate system. Throughout 
this work the 7 and 4 axis of the fixed coordinate axis will 

be directed east and north, respectively. Thus, the angle (/ 

is the angle between the wind direction and the east. 
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Figure 1. 
Coordinate Svstems Used in Derivation 


a af partially with respect to time gives 


. 2 ] 
=) Da a, a ae: +e c a a 
Sif il = cial MCS ca lao deal Sr TAG et Se 


or \? 
| ~<f sf zinl’ fez Zt ay | 
ay (af j: Be ied oe etsy szjeay~ by as Ly 6 f 
This may be written = eee, 
SVs gD (af) + 4 82 Be 
— 37 ie /* Cees Gd OT 
The expression on the right may then be substituted for 
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-2(7) in (2.3) and with the direct substitution for -/” 
in the first term on the right side results in 
ie : . ic] 
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Since the normal component of (2.1) may be written 


(FthVIV e422 2 0 





it is evident that the term inside braces in (2.4) is zero. 
Also, the left-hand side of (2.4) may be written in component 


form as follows: 


The tangential components of the left and right sides of 
(2.4) are thus identical, and the right side adds nothing which 
will assist in evaluating the isallohyptic gradient. while 
there may be some method of evaluating the tangential compon- 
ent of (2.3) without a direct knowledge of the time rate of 
change of the acceleration of the wind, it is not obvious. In 
view of the fact that acceleration is usually ignored in atmos- 
pheric motion, it might appear that the term =( ie is negli- 
gible. If this is ignored, the tangential component of (2.3) 
is 7 . 
Cases): Oberst (2.5) 
A qualitative interpretation of this equation with reference 
to the upper air chart (Chapter III) is not too difficult, 
and an err al check of it is discussed in Chapter IV. 

The normal component of (2.4) is 

eel ae it ey 1 

aa A llama oS ca A /-' 6S | (2.6) 

Except for the term involving a, all terms on the right- 
hand side of this equation can be svattated from a synoptic 
chart. In order to explore the possibility of substituting 
another expression for »h , Blaton's equation relating tra- 


jJectory curvature to streamline curvature is recalled. 
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It may be written i | ee : 


From this, taking the derivative with resvect to time, 
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Carrying out the indicated operation gives 


“is J “fp jf we? nd vi 
ae eee Sy eC e = Ma oe 
ofS a hpoy tia OW oy OP Se 48h corse 
' ’ H Ld a 
This may be written | » ae 
ie ose et Oy 
Poca) Ce? oe ory yt ye 


With this substitution, and after some rearrangement (2.7) 


becomes | 
Se eR, vet a+ vet sf -4ae 
ee fi — ie! a < aa 7 s ae + } pa : _ 
/ mer \, fee ~' ge \ ti ee yi Baa = Y d 
Pe 


But . it : ae. 


- is : eae 
LOY s ! oy f i = nwo — $ 
ee an =n +t |) nam \ - A a 7 \ } va 4 
cg ot / ae C ‘ y 
ana ali a ‘ aoe al 

Seale 6 Le : a Af eo al ry if s 

“, / Cv a on 
“id oe oe ee ee 
Cy { noe 


: 2a ay. poy! aF oY a4} 
Qh al fev ja yo BE + gt SS 4 Vee S 
Vx = Var (51 - Qn ot ary a5 TV Ss oe | (2.8) 


Substituting this back in (2.6) yields 
ry IC 3 oz oy wo 
74 Sn (SE ) = ey aky ar 935 t 82 (f- &KV) 


7 -O i re of | 4 ad 
VAQN) VEER ak EB 
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The tangential component of the isallohyotic gradient may be 


taken in its complete form from (2.3). It is 
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praGrae 
These two equations, (2.9) and (2.10), are the equations 
for the isallohyptic gradient normal to and along the wind 
vector, respectively. Equation (2.9) has been put in this form 
to make numerical evaluation of the equation from the weather 
map easier. There is some question whether this form is better 
than (2.6). It is to be noted that as written the units of 


"2 
~ 


Equations (2.9) and (2.10) are m sec ”, 





CHAPTER lil 


A QUALITATIVE INTERPRETATION OF THE 
ISALLOEYPTIC GRADIENT EQUATIONS 


Although (2.9) represents the form in which the normal 
component is evaluated empirically (Chapter IV), it is easier, 
perhaps, to gain an insight into the meaning of the equation 
from (2.6). 

The term qa (et), if positive at a given point on an 
isobaric surface, indicates that, relatively, the isobaric 
surface is rising more rapidly to the left of the wind direc- 
tion than to the right. Since the isobaric surface is lower 
to the left of the wind vector in virtually all instances 
above the friction layer, a positive value of this term indi- 
cates a relaxing isohyptic gradient. Similarly, a negative 
value of thie term can be associated with a tightening iso- 
hyptic gradient. 

The tangential term 95-(2F) of (2.10) is subject to a 
Similar interpretation. If it is positive, the isobaric sur- 
face is rising more rapidly (or falling less rapidly) ahead 
of the wind vector than behind it. Conversely, if the term 
is negative, the isobaric surface is falling ahead of the wind 
vector and rising behind the wind vector relative to the point 
where the term 1s being evaluated. 

In order to evaluate (2.10), it is necessary to evaluate 
separately the terms (th V2 and ¢(1¥) . The second tern, 


a /dV} 
S77 / , de difficult to evaluate, since it 1s the local rate 
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of change of acceleration. i1t is quite small, probably never 
exceeding 10°’ m sec™>. On the other hand, the first term 
(F+kv)V2e is also quite small. aA method for evaluating the 
latter term numerically is discussed in Chapter IV; however, 
the sign of this term may be determined by inspection in par- 
ticular cases. 

The factors / ana F¢ iV) are virtually always positive 
in the Northern Hemisphere. Thus, the sign of the first term 
on the right-hand side of (2.10) is determined completely by 
the sign of the factor oY . ina system of progressive open 
troughs and ridges such as that frequently found on upper- 
level charte, St is positive at the ridge line and negative 
at the trough line. This indicates that on the basis of this 
term alone, the tangential component of the isallohyptic gradi- 
ent points in the direction of the wind at the trough line 
and opposite to the wind at the ridge line, provided the trough 
and ridge are progressive. 

In onde for the normal component to be positive, it is 
necessary that the sum of the three terms on the right-hand 
side of (2.6) be greater than zero. None of these terms are 
inherently either positive or negative, and each must be exa- 
mined carefully in a given situation to determine its sign. 
These terms will be taken up singly in the following paragraphs 
to try to determine approximate magnitudes and sign. 

Let A= (fra: V4 3 a Se be an abbreviated notation for 
the first term of (2.6). This term consists of three factors, 


of which 4 is the only one which is always positive. The 


10 





factor £2 is positive if the wind is blowing across ischypses 


toward higher values. The factor $ is always positive in 

the Northern Hemisphere. The wind velocity VY is always 
positive. The curvature /{f may be positive (cyclonic curvature) 
or negative (anticyclonic curvature). The term es may be 
positive (counterclockwise turning) or negative (clockwise 
turning). The Coriolis parameter f is of the same order of 
magnitude as 2ik\// , though in most cases it is numerically 
greater. The term gy is usually not of the same order of 
magnitude as 7, although there are occasional exceptions. 

From this one can conclude that the factor (f p2KV 42) 
will usually be positive. An exception occurs when great 
anticyclonic curvature is associated with strong winds and 
the wind is locally turning clockwise with time. No such case 
has been found in the numerical evaluation of this factor. 

Thus, the sign of A is usually determined by the cross-— 
contour flow*of the wind. This term may vary in magnitude 
from 0 to 107° m sec™>. 

it is interesting to see what the contribution of this 
term is under various circumstances. Cross-—-contour flow to- 
ward higher values of 2 gives a positive contribution to 
2, (St) and indicates a weakening isohyptic gradient. Cross— 
contour flow toward low values of 7 indicates an increasing 
isohyptic gradient. The reverse of this would occur if 


(F 42h V3) were negative, which appears unlikely. 


ni! 





BA Pe osc 
The second term, B=V& (ftarv, ,» involves the factor 


(F724 V, . The variables appearing in this factor have just 
been discussed. Except for the case of abnormal anticyclonic 
flow (V>2V4 ) , this factor ia positive. The velocity,V , 

is positive alwavs. Therefore, the sign of this entire term 
is the same as the sign of the factor 2% in nearly all cases. 

The sign of 2 ids positive if the wind increases along 
the streamline (i.e. weak winds upstream with strong winds 
downstream) and negative if the wind is decreasing along the 
streamline. The magnitude of the term is a function of velo- 
city and rate of change of velocity along the trajectory in 
an obvious manner. It also is a function of curvature -- 
atrong cyclonic curvature being associated with greater magni- 
tudes, and strong anticyclonic curvature being associated with 
lower magnitudes. 

The term B alone then gives a vositive contribution in 
almost every case when weak winds are moving into an area of 
strong winds. This is the entrance region of Scherhag 9] : 
When such an entrance region is associated with cyclonic cur- 
vature, the contribution of B is positive and may be relatively 
large. An entrance region associated with anticyclonic curva- 
ture will almost always give a small positive contribution. 
This difference between the cyclonic and anticyclonic cases 
is again due to the factor (F+2KV) of 8 , discussed above. 

When strong winds are moving into a region of weak winds, 
the term § is nearly always negative. This is the delta 


region of Sherhag. A delta region associated with cyclonic 


iz 





curvature results in relatively large negative contributions 
to $i( #7) . The contribution of the delta region associated 
with anticyclonic curvature is negative, but relatively small. 
In magnitude this term also may vary from 0 to 107° m sec >, 


though it is doubtful that it approaches 107° 


very frequently. 
: i = 74d 
The third term on the right-hand eide of (2.6) is C=-/'Sh. 
Thus, if the curvature is increasing locally with time, term C 
gives a negative contribution. This would occur in a region 
into which a trough is moving. Conversely, a ridge moving 


fon 


into a region would give a positive contribution to a, (SF). 

Tne term C varies in magnitude from 0 to 107! m sec™?, 
and appears (See Chapter V) to be generally about one to two 
orders of magnitude smaller than term A . The range of mag- 
nitudes given for all three terms must be considered approxi- 
mate only. 

Several model-tvoe examples are now presented in which, 
Pacept ag noted, only the contributions due to terms A ana § 
are considered. These examples are illustrated in Figure 2. 

The first example depicted by Figure 2(a) is a situation 
in which both terms 4 and £& are positive. It should be 
noted that the winds, while flowing across contours toward 
higher values, are stronger downstream than upstream. There- 
fore, the isallohyptic gradient associated with this must be 
negative (i.e. a (3?) is positive), and the isohyptic gra- 


dient is relaxing. 


13 





The curvature of the trajectories does not affect fhe 
Sign of the isallohyptic gradient in this example unless the 
trajectory curvature is negative (anticyclonic) and exceeds 
the critical curvature for gradient flow. This case is dis- 
cussed below in connection with Figure 2(c). Trajectory cur- 
vature does change the magnitude of the isallohyptic gradient 
however, if all other things are held constant ~- cyclonic flow 
resulting in isallohyptic gradients of greater magnitude than 
anticyclonic flow. 

Figure 2(b) depicts a situation in which terms A ana DB 
are both negative. The chief features to be noted are the 
flow of air across the contours toward the low and the decrease 
of wind speeds downstream. The isallohyptic gradient asso- 
clated with this is positive and indicates that the isohyptic 
gradient is increasing with time. The influence of curvature 
is precisely the same in this case as for the first case de- 
scribed aboye. 

Figure 2(c) illustrates a condition often round at the 
crest of a ridge on upper-level charts. The curvature of 
the contours, given by Ro’, isa negative, but of greater 
magnitude than the maximum curvature, given by ee for 
gradient flow with the existing contour gradient. For this 
reason, air which is moving over the crest of the ridge will 
be unable to make the sharp turn required to remain in the 
contour channel, and so will cross contours toward lower iso- 


baric heights downstream of the ridge. If the trajectory of 


14 





~~ Bee a 
™ 
~ ae i“ 
~ a a 
Ba i a es, >» —_— =y 


Contours 


—_——-—-? Streamlines 


~ > Pra jectory 





Figure 2. Schematic Illustration of Synoptic Situations for 
Which Sign of Isallohyntic Gradient May Be Determined by 
Inspection. 
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such an overshooting air particle is a trajectory of maximum 
curvature, then the factor Ftibny; La =2ero. The proof or 
this statement 1s given below. 

Bjerknes / 1 | shows that for anticyclonic gradient flow 


« \/ _ 
lon = > . Va 


Now O37 4- ’ 


/ 
therefore, Rh, pe 
lway 


Hence #42 ieee VieaO 
This condition probably is met occasionally in the situa- 


tion of Figure 2(c). Then A=, and (2.6) becomes 


a (92) oz Ze 
9 Sul S¢) = a§ ; =) = ) (3.1) 


In this case it is not immediately obvious that the 
term | go is small compared with the other terms of tm 
equation, a8 has heretofore been assumed, nor are the two 
terms on the right side likely to be of like sign. Indeed, 
they are likely to give contributions of opposite sign. The 
factor oa will be positive, and ge will be negative, re- 
sulting in a negative contribution by the first term. The 
factor os will determine the sign of the second term. It 
will generally be negative, resulting in a positive contri- 
bution by the second term. The sign of the isallohyptioc 
gradient, then, is determined by the sign of the term of 


greater magnitude. 
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This 1s a synoptic situation which meteorologists usually 
agree is accompanied by greater rises of height to the right 
of the air stream than to the left. These rises take the 
form of a propagation of the crest of the ridge downstream. 
Thus, it appears that the negative, or first term on the right 
of (3.1) is the dominant one. In this example we may have one 
of the infreaquent cases in which oe is of the same order of 


magnitude as f. 


= 


1% 





CHAPTER LV 


NUMERICAL EVALUATION OF THE 
ISALLOHYFTIC GRADIENT EQUATIONS 


Obtaining general solutions of (2.9) and (2.10) is not, 
nBe general, feasible. It is possible, however, to obtain 
numerical solutions in particular cases, subject to certain 
simplifying assumptions. 

Equation (2.9) may be written in terms of finite differ- 





ences a a trajectory, as follows: ve 
Sm, es wy - 5) 
a Fr (22) = fo 2K, V)o (2274) pV ye (at It, V 
Avs AY, AAYs ahs (4.1) 
7 a (Zee 14) (FE) (- AW 


The subscripts refer to points located along a series of 


trajectories as shown in Figure 3. 
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rigure 3. Trajectories Showing the Grid Arrangement Required 
to Obtain a Numerical Solution 
of Equations (4.4) and (4.5). 
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Muoh of the success of evaluation depeiae on the accur- 
acy of the trajectories. Because of the fact that adequate 
upper-air charts are available only at 12-hour intervals, 
trajectories drawn by the method of successive approximation 
are of dubious value. Also, these require that the charts 
be available at both ends of the 12-hour time intérval. This 
makes the result useless as a prognostic tool. It was de- 
cided to make use of constant absolute vorticity trajectories 
(CAVT's) for the numerical work. 

The method of constructing CAVT's from inflection points 
ia quite generally known. The assumptions involved (non- 
divergence, constant lateral shear, frictionless flow, and no 
vertical motion) are quite restrictive. Equally restrictive 
for the vurposes of this study is the requirement that the 
CAVT be constructed from an inflection point. Martin [5 | 
ina study.which was carried out concurrently with this study 
pointed out,how CAVT's might be constructed from any point. 

A modi ficerton of his method was used in this study and is 
described in Appendix 1. 

The use of a CAVT as the trajectory simplified the eval- 
uation of (4.1) in a number of ways. The are length AS is 
equal to the product of V and AC . The curvature,i,, of 
the trajectory at point 2 (Figure 3) may be determined from 
the initial conditions and the latitude of point 2. Also, 
since the wind direction is known initially, it is a relatively 


simple matter to determine the turning of the wind at a given 
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point by measuring the angle between the initial wind direction 
and a tangent to the CAVT at that point. 
Equation (4.1) adapnted to the use of 12-hour CAVT's be- 


comes 


9 £5 (i= Je (tytn rg Ze2d , ous A (f+2 2i,V- 5% ) 
- ¥ (44 - 4m) _v*(Se)( tet a 


(4.2) 


Equation (2.10) written in similar form for numerical 


evaluation ia 
NPR. 
945(@) = (Fthavy - £155 (4.3) 


The second term on the right in (4.3) cannot be determined, 
but all other terms in (4.2) and (4.3) may be determined from 
a set of trajectories like those shown in Figure 3. Equations 
(4.2) and (4.3) with unit conversions, which permit data 
taken directly from an analyzed upper air chart to be used 


without change, become 


eu (Az) i, BQH sings t C053 hi: y)eotee ey 
A 


| (4.4) 
$0,594 (Vo V, (0524 20itg FOC IgV - SES 
9 208 16 Abe BAY) ERK IC™ “anf tet) 
and 
vAN A oe ee y 
aoe A 5a aay ee Jf CE AAI V AY 
A(A2) = 6.196 (0524 cond at £°) VMs waa 


The difference AZ found on the left side of these 
equations is the 12-hour change of height of the isobaric 


surface. The units in which the various parameters are ex- 
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preased are as follows: \/ (knots), H (reciprocal of degrees 
of latitude), Z4 (feet), Y and AY (degrees), 4” and JS 
(degrees of latitude). The isallohyptic ascendant is expressed 
in terms of feet per degree of latitude per 12 hours. All 
computed and observed values are of the isallohyptic ascendant. 

A series of 44 computations were made using CAVT's on the 
500 mb. level charts for 03007 and 12002 , 12 March 1954 and 
03002 , 13 March 1954. The only criterion used in selecting 
the points from which the CAVT's were constructed was that 
initial conditions be as clear-cut as possible. Approximately 
One-third of the trajectories were started from troughs. One- 
third were from inflection points on the west side of a trough, 
and the other one-third were from inflection points east of 
a trough. The fourth term of the right side of (4.4) proved 
to be so small that it was ignored. 

Another series of 21 computations were made from CAVT's 
on 500 mb. level charts selected from the period 9 through 18 
April 1954. The criterion used in selecting the CAVT starting 
point was that it be on an isotach maximum. This made the 
second term on the right in (4.4) negative in every case. 

All trajectories of both series have starting points within 
the United States. 

The computed normal components of the isallohyptic ascen- 
dant are shown for the series of 44, with the observed values, 
in Table 1. The coefficient of correlation between computed 
and observed values ig +0.65, indicating that approximately 


43% of the observed variation is accounted for by the computed 
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values. For 44 values the correlation coefficient, +0.65, is 
highly significant at the 99% level of belief, according to 
Hoel [ 4]. 

The regression equation (Figure 4 shows the scatter dia- 
gram and regression line of observed values ¥ on computed 


values XX.) is 


Y= B67 F080 4 | 

In this { and ¥ are expressed in units of feet per degree 
latitude per 12 hours. This equation indicates that in this 
series, the computed values were, for the most part, of the 
right sign, but too great in magnitude. The dominant term | 
in (4.4) was the cross—-contour term in most of the computa- 
tions in this series. An attempt was made to reduce this 
term empirically by multiplying it by a constant factor. This 
did not improve the correlation, however. | 

Table 1 also lists the computed and observed tangential 
components from this series of 44 trajectories. The correla- 
tion coefficient between these data is +t0.59, indicating that 
approximately 35% of the variance of the observed values is 
accounted for by the computed values. The correlation coef- 
ficient in this instance also is highly significant at the 
99% level of belief. 

The regression equation of observed values on computed 
values is 


Pea 5 3) eens 


(See Figure 5 for the scatter diagram and plot of this equation.) 
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TABLE 1 
Components of Isallohyptic Ascendant Associated 
with CAVT's Constructed on a 500 Mb. Chart 


Isallohyptic Ascendant 
(feet per deg lat per 12 hrs) 








Date 


Normal Component ‘Tangential Component Remarks 


Computed | Observed 





me a 
~ 23 


26 | southeast to- 


ce trough 


{ 


rrom inflec- 


tion point 







\ 


From trough 


to : northeast to- 
15002 | ward inflection 
12 | anne 
195k | | ; | 








ho 0 ) +120 | +33 | 

| + 14 | pe 0 From inflec- 

: =O pare al 8, + 28 ae point 

| — 21 + 43 + 26 northeast ES- 

| -12 || + 32 ward ridge 
15002 =. We Oo | 
12 Mar. +91 — 120 -— 37 iFrom inflec- 
to 0300 t 67 —~ 75 = 17 tion point 
13 Mar. Sie ie ~ 33 southeast to- 
1954 + 50 ~ 14 ~ 12 ward trough 





TABLE 1 
(Continued) 


Isallohyptic Ascendant | 
(feet per deg lat per 12 hrs) | 















em Normal Component || Tangential Component Remarks 
Time Computed | Observed Computed | Observed 
+85 | + 140 - ho ~ 46 Zz 
~ 200 - 91 + 140 — 38 r= trough 
=O -— 45 + 15 —~ 91 northeast to- 
15002 + 34 — 5 — 39 ~ 4s ward inflection 
m2 var, | + 27 0 — 36 ~ 56 point 
to 03008 — 76 - 62 +38 | +45 
Pa Mar. | -— 130 - 43 te | + 23 From inflec- 
1954 — 65 = mEO + 24 + 29 tion point 
| 4+- 34 a. + 27 abit. to- 
31 ward li 
a ma From inflec- 
+64 | + 80 ition point 
+ o57 + 62 southeast to- 





ward trough 


From trough 
northeast to- 


iward inflection 


point 

From inflec-— 
tion point 
northeast to- 


ward ridge 
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Figure 4. Scatter Diagram and Regression Line of Computed 
against Observed Normal Isallohyptic Ascendant. 
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Figure 5. Scatter Diagram and Regression Line of Computed 
against Observed Tangential Isallohyptic Ascendant. 
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TABLE 2 


Components of Isallohyptic Ascendant Associated 
with CAVT's Constructed from Isotach 
Maxima on a 500 Mb. Chart 


Isallohyptic Ascendant 
(feet per deg lat per 12 hrs) 











| 
Date | 
and | Normal Component Tangential Component, Remarks* 
Time ae me Observed Site ag Observed 
1954 a 3 en | 7 
3/9-03 | Sie ~ 5 + 4 | Intl. to vers 
3/9-15 | - 170 0 +2h 0 Near infl. to 
| : | | lie. 
3/10-03| — 3 Hom | 0 +29 |Rag. to tr. 
| 
3/10-15 | — 33 — Lo + 34 +17 | inf. tO tre 
3/11-03; ~- 210 + 38 + 160 +43 ler, to. ings 
| 
B/11-15 | + 28 + 14 | — 9 | + 14 Rag. CO. ce 
ali-15| — 20g | +55 || +23 +40 |Tr. arnd. top 
| | of low 
: 
o2-03! -“g@i8O | -— 100 0 + 46 Tr. arnd. top 
: | of low 
| - 
3/12-15) =F 0 0 + 25 Tr. apna top 
| I of low 
3/13-03} -160 | +16 || +200 + 9 
| | 
4/10-15; -130 | 2 Oe + 43 — 33 Tr. to. 2g 
| : 
4/11-03| ~370 | Oo | + 230 0 Tr. to. Dini 
4/13-03; + 33 © - 8 | me ike | + 46 Infl. to. Ge 
b/1h=15 | = Oo § + 20 + 30 — 17 itr. to. rdg. 
4/15-03| -42 4 — 23 — 8 — 14 'Infl. to. gen 
4/15-15; +190 — — 17 — 150 — 54 Near infl. to 
. tr. 
4/16-03] + 150 +150 — 68 - 19 Infl. to tm 





TABLE 2 
(Continued) 


Isallohyptic Ascendant 
(feet per deg lat per 12 hrs) 





Date 
and Normal Component Tangential Component Remarks* 
Time Computed | Observed Computed jObserved 
| 
whe 15) +180 —- 56 — 210 —12 ss to tre 
4/17-03| ~— 180 +140 +110 +12 ‘Near tr. to 
| | rdg. 
4/17-15) + 95 - 95 + 115 ey Infl. acr. wk. 
rag. andes 
4/18-03{} + 120 — 40 0 0 ‘ag tr. sen. 
wk. rdg. and 
sin tr. 


*Abbreviations used in this column are: Infl. inflection point, 
tr. trough, rdg. ridge, wk. weak, str. strong, arnd. around, 


acr. across, to. toward. 


* 


“O 
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As in the case of the normal component, this indicates that 
comouted values were predominately of the right sign, but of 
too great magnitude. 

Of the set of 21 trajectories which were constructed from 
{sotach maxima, less than half resulted in computed normal 
components of the right sign. Also, as may be observed from 
Table 2, the magnitude of the computed values were almost 
invariably larger than the magnitude of corresponding observed 
values. 

The correlation coefficient between observed and computed 
values is -0.09, which is not significant at any acceptable 
fevel of belief. 

The computed tangential components resulting from this 
set of 21 trajectories were somewhat better than the computed 
normal components. These (Table 2) were found to have a cor- 
relation coefficient of +0.41 with observed values. This is 
significant-at the 93% level of belief. Since it indicates 
that the aut values account for only 17% of the variation 
occurring in the observed data, it cannot be considered a good 
forecast tool. In this case, as in the others, the magnitude 
of the computed values was too large. No regression equation 
was developed for either of the components computed from this 
series of trajectories. 

Because of the large contribution made by the cross- 
contour term in most of these computations, the CAVT's were 


examined to try to determine what might be done to the trajec- 
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tory to reduce this contribution. Usually, though not always, 
a reduction of the initial wind speed would have resulted in 

a decrease in the contribution of the first term of (4.4), 
without changing appreciably the contribution due to the other 
terms. This might have improved the over-all results; however, 
more study is required before a definite statement can be made 
in this regard. 

In using the difference equations, (4.4) and (4.5) instead 
of the differential equations, (2.9) and (2.10), it is assumed 
that the differences are small enough so that the difference 
equations closely approximate the integrated forms of the dif- 
ferential equations. This assumption is almost surely not 
valid when the time interval is 12 hours. Even with a perfect 
12-hour trajectory, the results would be somewhat uncertain 
because of this. 

Another source of error is found in the crudeness with 
which the bagic data must be obtained. Curvature plays an im- 
portant sented this problem in several ways; yet it is diffi- 
cult to determine the curvature of the contours on an upper 
alr chart with any degree of accuracy. The curvature of the 
air trajectory is even less attainable. Other parameters are 
equally difficult to evaluate. This is especially so in con- 
nection with short waves for which the measurement of the 
initial parameters becomes increasingly critical. This pre- 
sumably accounts for some of the lack of agreement of Table 2. 


Furthermore, very severe verification criteria were imposed: 
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The computed normal isallohyptic component was verified along 

& curve roughly orthogonal to the center point of the various 
CAVT's. It is quite possible that only a slight displacement 
of the forecasted isallohyptic centers located by the computa- 
tional procedure would lead to the variation between the ob- 
served and computed isallohyptic gradients. The best test of 
the 12-hour forecast procedure would be to integrate the 
isallohyptic tendency equations to determine the complete fore- 
cast isallohyptic pattern, and then verify the location of rise 
and fall centers against the observed pattern. The forecasted 
isallohyptic gradients correlate sufficiently well (see Chapter 
V) with the observed gradients, in spite of the severe verifi- 
cation just mentioned. 

In certain of the cases presented in Table 2, the CAVT 
frequently was constructed from a minor trough or ridge or an 
inflection point between minor features. In such cases, the 
CAVT frequently crossed a synoptic feature such as a ridge (or 
trough) , continuing right on part way across the next trough 
(or ridge) downstream. Trajectories of this type may be re- 
latively unrepresentative. The best test in such cases, also, 
would be that of verification of the forecast isallohyptio 


patterns. 
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CHAPTER V 
CONCLUSIONS 


in Chapter III it was shown that intensification of cer- 
tain synoptic models on an isobaric surface can be determined 
qualitatively, provided the required parameters are determin-— 
able. Of the parameters required, the most difficult to deter- 
mine are the cross-contour flow and the local time rate of 
turning of the wind. The cross-contour flow 1s by far the more 
important of these two parameters when using (2.6) or (2.9) 
with trajectories as drawn in this dissertation. 

O'Connor's hypothesis fé] seems to indicate that the most 
important term in (2.6) is the second term on the right-hand 
side, since his discussion points out, primarily, those features 
associated with term 8. This follows since O'Connor's point 
of view is that an isotach maximum or minimum will remain 
within its &ppropriate contour channel at each stage of its 
movement. Thus, no cross-contour flow occurs in this scheme. 
Furthermore, in many cases, O'Connor's method seems to give 
acceptable qualitative results. 

With the tvpe of trajectories used in this work, the 
verification was also quite good, at least in the data of 
Table 1. As discussed in the next paragraph, the cross-contour 
term A is, in general, much larger than the longitudinal wind 
difference, term B , and the former gives an over-estimate 


(with correct sign) of the normal component of isallohyptic 
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gradient. Thus, the particular objective technique employed 
herein, while demonstrating the importance of term Ai , appears 
to over-estimate it and under-estimate term © . Both of these 
shortcomings might well be overcome by extracting data from a 
forecast map halfway through the forecast period. This would 
have required application of an iterative forecast procedure 
involving time intervals of 1 - 2 hours at a step. However, 
the scope of the vresent work was to devise an objective fore- 
cast procedure using only the initial map. 

Some information concerning the relative importance of the 
various terms of (2.9) may be gained by further study of the 
more successful series of computations, the 44 which are sun- 
marized in Table 1. The average contribution of each of the 
terms (averaged without regard to sign) are: term A (cross- 
contour term) 73.4, term 8 (longitudinal wind difference 
term) 5.7, term © 2.2. The fourth term was ignored because 
it was smali.,compared with the other terms. The average ob- 
served uienate 34.8, and the average total computed value is 
78.1. Further, the coefficient of correlation between the 
values of the first term alone and the observed values is 
t0.55. This is significant at the 99% level of belief. The 
coefficient of correlation between the second term values and 
the observed values is -0.03. This is not significant at any 
acceptable level of belief. Also not significant is the 
-0.04 correlation coefficient between the computed values of 


terms A and B . 
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Although terms 6 and CG of (2.9) are small, there is 
some irprovement in the correlation between observed and com- 
puted values when they are included. The discussion at the 
beginning of this chapter seems to indicate that there must be 
a certain correspondence between the signs of the first and 
second terms. Thus, if the second term contribution can be 
used to indicate the sign of the entire computed value, even 
though its value is only a small fraction of the total value, 
it should have the same sign as the contribution from the first 
or largest term. This is not borne out by the data of this 
series. Of the 44 pairs of values obtained by matching cor- 
responding terms A ana B , there were 38 pairs having non- 
zero second term values. These 38 pairs consisted of 20 pairs 
having opposite sign and 18 pairs having the same sign. This 
lack of correspondence of signs is also indicated by the —0.04 
correlatiom coefficient between these values. 

16e) conelusion it appears that, while significant results 
have been are eA by using CAVT's drawn from selected points, 
some more reliable method of determining air trajectories is 
necessary to make this method a really useful prognostic tool. 
The evaluation of all terms of (2.9) and (2.10) is dependent 
On an accurate air trajectory, but it is in the evaluation of 
term A of (2.9) that trajectory errors are most keenly felt. 
Without actually constructing trajectories, however, some in- 
formation can be obtained by qualitative application of the 


principles which have been set forth. The examples discussed 
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in Chapter III represent synoptic patterns which are amenable 


to qualitative interpretation. 
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APPENDIX I 
& METHOD OF CONSTRUCTING CAVT'S FROM ANY POINT 


After making the usual assumptions (i.e. non-divergence, 
constant lateral shear, horizontal frictionless motion), the 


absolute vorticity equation may be written 
Vi Fame 24 ys 
7. 1 le (1) 
The subscripts indicate initial conditions. 
Subject to the usual additional assumption that V is 
constant, 4 may be determined at any latitude for a given set 


of initial conditions. Thus 
VA 


or 


a ee wa 
of \ pLr- (2017 PAE? q)+ Ko 


Substituting appropriate constants and conversion factors 


e 


to permit the use of VY in knots and 77, and F in degrees 


latitude changes this to 
a 
4 (an o-ant) + 
Radii may be computed for the desired latitudes (either 
higher or lower than the initial latitude, depending on the 
initial wind direction). Successive arcs of circles drawn 


end on end using appropriate radii result in the CAVT. A 


CAVT constructed by this means is not a curve of continuously 
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varying curvature. By increasing the number of arcs used in 
constructing the CAVT, it may, however, be made to approach 
as closely as is desired the "true" CAVT which has continuously 
varying curvature. Generally, the computation of radii at 
5 degree latitude intervals is adequate. 

Figure 6 illustrates the actual CAVT construction process. 
It is to be noted that the center of curvature of each succes- 
sive arc lies on a line drawn through the starting point of the 
arc and the center of curvature of the last arc. This insures 
continuity of the curve between arc segments. Table 3 lists 
the initial and computed data used in the CAVT conetrueimon 


shown in Figure 6. 


TABLE 3 


Initial and Computed Data Used in CAVT 
P Construction Shown in Figure 6 
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Figure 6. Constructing a CAVT from a Point Other Than an 
Int tection= rein. 
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